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ABSTRACT 

Inventory is essential for the efficient running of any business. Single location inventory system is considered by 
many researchers. This paper deals with two echelon inventory system with handling two products having joint ordering 
policy. The demand for the products follows independent poison distributions at retailer and distributor node. The items are 
supplied to the retailers from the distribution center (DC) administrated with exponential lead time having parameter (J. 
(>0). The joint probability disruption of the inventory levels of two products at retailer and the supplier are obtained in the 
steady state case. Various system performance measures are derived and the long run total expected inventory cost rate is 
calculated. Several instances of numerical examples, which provide insight into the behavior of the system, are presented. 

KEYWORDS: Continuous Review Inventory System, Two-Echelon, Joint Ordering Policy 

1. INTRODUCTION 

Inventory control system in two-echelon or multi echelon has been considered by many researchers. Most of the 
papers the authors concentrate on periodic review inventory system. A first quantitative model to find an optimal (s, S) 
inventory policies for single item at single location is considered by Veinott and Wagner[10]. Thangaraj and 
Ramanarayanan[9] develop an inventory model with two re-order levels and random lifetime. Gross and Harris [7] 
considered the inventory system with state dependent lead times. All the above papers deals only with single location 
inventory system. This paper deals with two echelon inventory system(retailer -lower echelon, distributor - higher 
echelon) handling two different products supplied by the same distributor. 

Anbazhagan and Arivarignan [1] have analyzed two commodity inventory system under various ordering policies. 
Yadavalli et. al., [11] have analyzed a model with joint ordering policy and varying order quantities. Yadavalli et. al., [12] 
have considered a two commodity substitutable inventory system with Poisson demands and arbitrarily distributed lead 
time. 

In a very recent paper, Anbazhagan et. al. [2] considered analysis of two commodity inventory system with 
compliment for bulk demand in which, one of the items for the major item, with random lead time but instantaneous 
replenishment for the gift item are considered. The lost sales for major item is also assumed when the items are out of 
stock. The above model is studied only at single location(Lower echelon). We extend the same in to multi-echelon 
structure (Supply Chain) with joint ordering policy. The rest of the paper is organized as follows. The model formulation is 
described in section 2, along with some important notations used in the paper. In section 3, steady state analysis are done: 


www.iaset.us 


editor @iaset.us 


44 


R. Karthikeyan & K. Krishnan 


Section 4 deals with the derivation of operating characteristics of the system. In section 5, the cost analysis for the 
operation. Section 6 provides Numerical examples and sensitivity analysis. 

2. MODEL 

2.1. The Problem Description 

The inventory control system considered in this paper is defined as follows. We assume that finished products are 
supplied from warehouse to distribution centre (DC) which adopts (0, M) replenishment policy then the product is supplied 
to retailer (R) who adopts (Sj, Q,) policy. The demands at retailer node follows independent Poisson distribution with rate 
A-i (i = 1, 2). The items are supplied to the retailers in packs of Q (Q = Q 1 +Q 2 ) where Q, (= Sj-S;) from the distribution center 
(DC) administrated with exponential lead time having parameter |X (>0). The direct demand at distributor node follows 
Poisson distribution with rate Ad- The replenishment of items in terms of pockets is made from WH to DC is instantaneous. 
Demands that occur during the stock out periods are assumed to be lost sales. In this model the maximum inventory level 
at retailer node Si is fixed and he reorder level is fixed as Sj for the i-th commodity and the ordering policy is to place order 
for Qi (= Sj-Sj) items (i = 1,2) when both the inventory levels are less than or equal to their respective reorder levels. The 
maximum inventory level at DC is M (M = 11 Q). The joint probability distribution for both commodities is obtained in 
steady state cases. 

The optimization criterion is to minimize the total cost rate incurred at all the location subject to the performance 
level constrains. According to the assumptions the on hand inventory levels at all the nodes follows a random process. 

2.2. Analysis 

Let Ij(t), (i = 1, 2,3) denote the on-hand inventory levels for commodity- 1, commodity-2 at retailer node and 
Distribution Centre (DC) respectively at time t+. From the assumptions on the input and output processes, I(t) = { Ij(t); t> 
0 } (i = 1, 2,3) is a Markov process with state space 


E = 


(i, k, m ) / i = Sj , (Sj — 1), ... , sj , (sj - 1), ...,2,1,0.,' 
< k =S 2 ,(S 2 -1), ,s 2 ,(s 2 -1), ... ,2,1,0. 

m = nQ,(n-l)Q, ... ,Q 


Since E is finite and all its states are recurrent non-null, I(t) = { Ij(t); t > 0 }is an irreducible Markov process with 
state space E and it is an ergodic process. Hence the limiting distribution exists and is independent of the initial state. The 
infinitesimal generator of this process 

R = (a(i,k,m : j,l,n)) (ikm) (j a n)eE 


can be obtained from the following arguments. 

• The arrival of a demand for commodity- 1 at retailer node makes a state transition in the Markov process from 
(i, k, m) to (i-1, k, m) with intensity of transition / H . 

• The arrival of a demand for commodity-2 at retailer node makes a state transition in the Markov process from 
(i, k, m) to (i, k-1 m) with intensity of transition A 2 . 
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• The arrival of a demand at distributor node makes a state transition in the Markov process from (i, k, m) to (i, k 
m-Q) with intensity of transition Ad- 

• Joint Replenishment of inventory for commodity- land commodity-2 at retailer node makes a state transition in the 
Markov process from (i, k, m) to (i+Qi, k+Q 2 , m-Q) with intensity of transition p .where Q = Q 1 +Q 2 

The infinitesimal generator R is given by 

B 0 ••• 0 0 

0 A B ••• 0 0 

0 0 A ••• 0 0 

R = 

0 0 0 ••• A B 

V B 0 0 ••• 0 A 

The entities of matrix R are given by 

A if p = q ;p = nQ,(n-l)Q, ... , 3Q, 2Q. 

B if p = q+Q ; P = nQ ,(n-l)Q, ... , 3Q, 2Q. 

B if p = q-(n-l)Q ; p = Q. 

0 otherwise 

The sub matrices of matrix R are given by 

A x if p = q ; p = S 2 ,(S 2 -1),(S 2 -2), ... ,(s 2 +l). 

A 2 if p = q+Q ; p = S 2 ,(S 2 — 1),(S 2 -2), , 1. 

[ A lpxq =<A3 if P = q ;P = S 2> (s 2 - 1 )»((s 2 -2), - I- 
A4 if p = q ; p = 0 
0 otherwise 

p if p = q + Q 2 ; q = s 2 , (s 2 -1), ... ,1,0 
0 otherwise 

The sub matrices of matrix A are given by 

ifp = q-l ;q = 1,2, ... ,S l 
-(A, + X 2 + X D ) if p = q-1 ; q = 1,2, ... , S : 

-(A, 2 +X D ) ifp = q;q = 0 

0 otherwise 
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J A 2 if p = q-l ; q = 1,2, ... , S : 

[ 0 otherwise 

\ ifp = q-l ;q= 1,2, ... ,S 1 

-(^ + \ 2 + + 1 1 ) ^ p = q _ i ; q = 1,2, ... , s : 

-(A 2 + A D + |l) if p = q ; q = 0 
0 otherwise 


\ if p = q-1 ;q= 1,2, ... ,S l 

-(Aj + A 2 + A d + q) if p = q-1 ; q = 1,2, ... , Sj 

-(|i + A D ) if p = q ; q = 0 


0 otherwise 


4.2.3 Transient Analysis 

Define the transient probability function 

P<i,k,m)(j, 1, n : t) = Pr {(1,(1), I 2 (t), I 3 (t)) = (j, 1, n) I (1,(0), I 2 (0), I 3 (0)) = (i, k, m)}. 

The transient matrix for t > 0 is of the form P(t) = tp n k mJ (j, 1, n : t)) (i k m)( j , n)eE satisfies the Kolmogorov- 
forward equation P (t) = P(t).R , where R is the infinitesimal generator of the process {I(t), t > 0 } . From the above 

equation, together with initial condition P(0) = I, the solution can be express of the formP(t) = P(0)e^ 1 =e Rt , 
where the matrix expansion in power series form is 


e 


Rt 


oo 


I+Z 

n=l 


rV 

n! 


4.2.4 Steady State Analysis 

The structure of the infinitesimal matrix R reveals that the state space E of the Markov process {I(t);t > 0} is 
finite and irreducible. Let the limiting probability distribution of the inventory level process be 

=lim Pr { (i 1 (t),i 2 (t),i D (t))=(i,k,m) } , 

t— »oo 


where m rk is the steady state probability that the system be in state (i,k,m), (Cinlar [5]). 

Let n = ( nQ n, (n-1)Q II, <n 2)Q II, 2Q II, Q II. ) denote the steady state probability distribution, 
where l<-) II = ( 1 , 1 1 ^ , . . . , II J' , II Jj j for j = 1, 2 ... n and k = 1 , 2 ,..., S for the system under consideration. Lor each 
(i,k,m), m nf can be obtained by solving the matrix equation II R = 0 

Since the state space is finite and R is irreducible, the stationary probability vector II for the generator R always 


www.iaset.us 


editor @iaset.us 


Analysis of two Echelon Inventory System with Joint Ordering Policy 


47 


exists and satisfies 

IIR=0 and 1 1 e= 1 

The vector II can be represented by 

_ ^j-j<Q1,Q2> j-j<2Q1,2Q2> j-j< 3Q1,3Q2> j-j<nlQl,n2Q2> 

Now the structure of R shows, the model under study is a finite birth death model in the Markovian environment. 
Hence we use the Gaver algorithm for computing the limiting probability vector. For the sake of completeness we provide 
the algorithm here. 

4.2.5 Algorithm 

1 . Determine recursively the matrix Dn, 0 < n < N by using 

DO = A 0 

Dn= A n +B n (— D n l_i) C, n = 1, 2,....K 

2. Slove the system 

n <N> D N =o 

3. Compute recursively the vector n <n> , n = N — 1, ,0 using 

n<n> = n <n+ 1 > B n+ 1 (-D“ 1 ),n = n-l, ,0 

4. Re -normalize the vector n, using 

rie = l. 

4.2.6. Operating Characteristics 

In this section, we derive some important system performance measures. 

(a) Mean Inventory Level 

Let IL R1 and IL R2 denote the expected inventory level in the steady state at retailer node for the commodity- 1 and 
commodity-2, and IL D denote the expected inventory level at distribution centre. They are defined as 


nQ 

S 2 Sl 


il ri = X 

* ^ ^ •j-j«<i,k,m.o»> 

(4.1) 

m=Q 

k=0 i=0 


n Q 

S 2 Si 


IL R 2= ^ 

* ^ ^ j,j-j<«i,k,m,o>» 

(4.2) 

m=Q 

k=0 i=0 



Impact Factor (JCC): 2.6305 


NAAS Rating 3.19 



48 


R. Karthikeyan & K. Krishnan 


nQ * S 2 S, 

il d= £ ££mn <<<1 ’ k -° >>> ;Q = Q 1+ Q 2 (4.3) 

m=Q k=0 i=0 


(b) Mean Reorder Rate 


Consider the reorder events p R for both commodities at retailer node and p D of at distribution centre. It is observe 
that p D event occur whenever the inventory level at DC node reaches 0 whereas the event p R occurs whenever the 
inventory level at retailer node drops to either (si,s 2 ) or (si,j),j < s 2 or (i,s 2 ),i < s R 

The mean reorder rate at retailer node and distribution centre are given by 


nQ f s 2 Sj 

Br= X X^i) n<<<<Si+1 ’ k,m>>>>+ X^2) n<<<i ’ S2+1 ’ m>>> 

m=QVk=0 i=0 


(4.4) 


» 2 * 1 


b d = EZ(P^ D )n <<<<lW>>>> 


(4.5) 


k=0 i=0 

(c) Mean Shortage Rate 

Shortage occurs only at retailer node and the mean shortage rate at retailer is denoted by a R which is given by 


n Q ^ f $2 S] ^ 

«R = X X^l) n <« °’ k ’ m >» +X ^2 ) n «<i,0,m>» 

m— Q Vk-0 k-0 ) 


(4.6) 


4.2.7 Cost Analysis 

In this section we analyze the cost structure for the proposed models by considering the minimization of the 
steady state total expected cost per time. 

The long run expected cost rate for the model is defined to be 

C(s,Q) — h RI IL R1 + h R2 IL R2 + h D IL D + k r B R +k D B D + g R ot R (4.7) 

where 

• h R1 , h R2 , and h D are the holding cost per unit of item of Commodity- 1, Commodity-2 at retailer node and holding 
cost per unit items (Q items, Q = Qi+Q 2 ) at distribution centre respectively per unit time. 

• ILri, IL R2 , and IL D are the inventory level of Commodity- 1, Commodity- 2 at retailer node and inventory level 
at distribution centre respectively per unit time. 

• K r , and K D are the Fixed ordering cost at retailer node distribution centre respectively. 

• B r and B d are the mean reorder rate at retailer node and distribution centre. 

• a R is the mean shortage rate at retail node. 
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• g R is the shortage cost per unit shortage at retailer node. 

Although, we have not proved analytically the convexity of the cost function TC(s Iy s 2 Q) ,our experience with 
considerable number of numerical examples indicates that TC(si,S 2 Q) for fixed Q to be convex in s. In some cases, it 
turned out to be an increasing function of s. Hence, we adopted the numerical search procedure to determine the optimal 
values s*, consequently, we obtain optimal n*. For large number of parameters, our calculation of TC(s h s 2 Q) revealed a 
convex structure for the same. 

5. NUMERICAL ILLUSTRATION 

In the section the problem of minimizing the long run total expected cost per unit time under the following cost 
structure is considered for discussion. The optimum values of the system parameters s t and s 2 are obtained and the sensitive 
analysis is also done for the system. 

The results we obtained in the steady state case may be illustrated through the following numerical example, 
S,=20, S 2 = 25, M= 200, ^=3,^=4^ A D = 2 , fj, =3 

K= 1.1, h D =1.2, k R =l.5,k d =1.3 g R =2.1, g D =2.2, C 0 =2.3 


The cost for different reorder level are given by 

Table 1: Total Expected Cost Rate as a Function Si,S 2 and Q 


Si 

2 

3 

4 

5* 

6 

7 

8 

s 2 

3 

4 

5 

6* 

7 

8 

9 

QMQ 1 .+ Q 2 ) 

40 

38 

36 

34 

32 

30 

28 

TC(s„s 2 0) 

115.7454 

105.8264 

90.5506 

87.90082* 

94.28986 

110.392 

123.6665 


For the inventory capacity S, the optimal reorder level s* and optimal cost TC (s 2y s2, Q) are indicated by the 
symbol *. The Convexity of the cost function is given in the graph. 



1.1. Sensitivity Analysis 

The effect of changes in Demand rate at retailer node for product 1 and 2. 

Table 2: Total Expected Cost Rate as a Function When Demand Increases 



>> 

II 

00 

O 

II 

<< 

II 

<< 

II 

<< 

'sO 

II 

<< 

>.!= 12 

45.073539 

46.419603 

47.214733 

47.629158 

47.785721 

^=14 

46.271834 

47.617898 

48.413028 

48.827452 

48.984016 
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Xi= 16 

47.136498 

48.482562 

49.277692 

49.692117 

49.848680 

II 

M 

00 

47.793130 

49.139194 

49.934324 

50.348749 

50.505312 

o 

<N 

II 

c< 

48.325884 

49.671948 

50.467078 

50.881502 

51.038066 


The graph of the demand rate variation is given below and it describes, if the demand rate increases then the total 
cost also increases. 



Figure 2: TC (S h S2, Q) for Different Demand Rates 


Table 3: Total Expected Cost Rate as a Function When Si and Si Increases 



Si=2 

Si=4 

S,=6 

S,=8 

S 1= 10 

Si=45 

111.998164 

118.761642 

124.799813 

130.664478 

136.843108 

Si=50 

118.699318 

124.359054 

129.299781 

133.996082 

138.773110 

Si=55 

123.090318 

127.562168 

131.434202 

135.090008 

138.765317 

Si=60 

125.576395 

129.128145 

132.243072 

135.192545 

138.137225 

Si=65 

126.984419 

130.065937 

132.851239 

135.511516 

138.145444 



Figure 3: TC (S h S2 , Q) For Different Si and Si Values 


From the graph it is identified that the total cost increases when the s and S increases. 

Table 4: Total Expected Cost Rate as a Function When Si and S 2 Increases 



S 2 =30 

S 2 =35 

O 

Tf 

II 

n 

S 2 =45 

S 2 =50 

S!=20 

189.8666 

179.2934 

165.6356 

152.2989 

141.219 

S 1= 25 

195.0113 

186.1366 

173.264 

159.6361 

147.8542 

S 1= 30 

199.5703 

192.4455 

180.7208 

167.1007 

154.6623 

Sj=35 

203.6145 

198.1655 

187.8799 

174.6251 

161.6628 

S 1= 40 

207.2345 

203.2862 

194.623 

182.1144 

168.8466 
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Figure 4: TC (S h S2 : Q) for Different Si and S 2 Values 

From the graph it is identified that the total cost decrease when S 2 increases and increases when Si increases 


Table 5: Total Expected Cost Rate When H r and H d Increases 



h D =0.04 

h D =0.08 

h D =0.12 

h D =0.16 

h D =0.20 

h R =0.002 

43.1924942 

43.87908 

44.5655 

45.25206 

45.93862 

h R =0.004 

43.2037231 

43.89028 

44.5767 

45.26326 

45.94982 

h R =0.006 

43.2149519 

43.90148 

44.58804 

45.27446 

45.96102 

h R =0.008 

43.2261794 

43.91268 

44.59924 

45.28566 

45.97222 

h R =0.010 

43.2374082 

43.92388 

44.61044 

45.297 

45.98342 



Figure 5: TC (S ,, S2 i Q) for Different H r and H d Values 

As is to be expected the graph shows that the total cost increases when h R and h D increases 


Table 6: Total Expected Cost Rate When G r and K d Increases 



gR-0.2 

gR-0.4 

g R =0.6 

g R -0.8 

g R =l-0 

FT 

O 

II 

o 

32.346268 

33.06786 

33.78945 

34.51104 

35.23264 

k D =15 

34.204663 

34.926255 

35.64785 

36.36944 

37.09103 

k D =20 

36.063058 

36.78465 

37.50624 

38.22783 

38.94943 

k D =25 

37.921453 

38.643045 

39.36464 

40.08623 

40.80782 

k D =30 

39.779848 

40.50144 

41.22303 

41.94462 

42.66622 
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Figure 6: TC (Si, 52 Q) for Different G R and K D Values 

As is to be expected the graph shows that the total cost increases when g R and k D increases. 

6. CONCLUSIONS 


This paper deals with a two echelon Inventory system with two The demand at retailer node follows independent 
Poisson with rate ^ and X 2 . The Joint replenishment ordering policy is applied at retailer node. The direct demand at 
distributor is assumed to be Poisson with rate X D The structure of the chain allows vertical movement of goods from to 
supplier to Retailer. If there is no stock for product at retailer the demand is refused and it is treated as lost sale. The model 
is analyzed within the framework of Markov processes. Joint probability distribution of inventory levels at DC and Retailer 
for both products are computed in the steady state. Various system performance measures are derived and the long-run 
expected cost is calculated. By assuming a suitable cost structure on the inventory system, we have presented extensive 
numerical illustrations to show the effect of change of values on the total expected cost rate. It would be interesting to 
analyze the problem discussed in this paper by relaxing the assumption of exponentially distributed lead-times to a class of 
arbitrarily distributed lead-times using techniques from renewal theory and semi-regenerative processes. Once this is done, 
the general model can be used to generate various special eases. 
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